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CSJ ■ ABSTRACT 

J> ' The population of high energy and very high energy gamma-ray sources, detected with 

ON , EGRET and the new generation of ground-based Cherenkov telescopes, conforms a reduced 

but physically important sample. Most of these sources are extragalactic (e.g., blazars), while 
among the galactic ones there are pulsars and supernova remnants. The microquasar LS 5039, 
previously proposed to be associated with an EGRET source by Paredes et al. (2000), has 
recently been detected at TeV energies, confirming that microquasars should be regarded as a 
\ class of high energy gamma-ray sources. To model and understand how the energetic photons 

— . i are produced and escape from LS 5039 it is crucial to unveil the nature of the compact object, 

*-£h ■ which remains unknown. Here we present new intermediate-dispersion spectroscopy of this 

source which, combined with values reported in the literature, provides an orbital period of 
Q ■ Porb = 3.90603 ± 0.00017 d, a mass function f(M) = 0.0053 ± 0.0009 M Q , and an eccen- 

| tricity e = 0.35 ± 0.04. Atmosphere model fitting to the spectrum of the optical companion, 

r/3 , together with our new distance estimate of d = 2.5 ± 0.1 kpc, yields Ro — 9.3+q'q Rq, 

log(L /L©) = 5.26 ± 0.06, and M = 22.9:^:9 M o- These, combined with our dynamical 
. solution and the assumption of pseudo-synchronization, yield an inclination i = 24.9 ± 2.8° 

■ and a compact object mass Mx — 3.7^*' q Mq. This is above neutron star masses for most of 

rN | the standard equations of state and, therefore, we propose that the compact object in LS 5039 

• is a black hole. We finally discuss about the implications of our orbital solution and new 

parameters of the binary system on the CNO products, the accretion/ejection energetic bal- 
ance, the supernova explosion scenario, and the behaviour of the very high energy gamma-ray 
emission with the new orbital period. 

Key words: stars: accretion, accretion discs - binaries: close - stars: individual: LS 5039 - 
X-rays: binaries - X-rays: individual: RX J1826.2— 1450. 
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1 INTRODUCTION 

The third EGRET catalog contains 27 1 sources detected at energies 
above 100 MeV jHartman et ailll999t) . Apart from extragalactic 
sources at high galactic latitudes and some galactic pulsars and su- 
pernova remnants, the majority of these sources, ~ 168 or ~62 per 
cent, still remains unidentified. Among them, there are 72 sources 
located at low galactic latitudes, having |&|<10°, which represents 
around 45 per cent of the unidentified sources. Therefore, several 
of these objects are presumably of galactic nature. The existence 
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of similar properties between some of these sources, indicate that 
there are at least three different groups of galactic popul ations: a 
group of young stellar objects and star-forming regions jRomerol 
2001), the sources forming a halo around the galactic cent er, and 
finally a gr oup of sources corre lated with the Goul d Belt iGrenieJ 
l200d) (see iRomero et ai]l2004l and iGreniejEool for recent up- 
dates). The identification of their nature is of prime importance in 
high energy astrophysics (see lCheng & Romerdl2005l for a recent 
compilation). 

On the other hand, microquasars are galactic X-ray bina- 
ries with relativistic collimated jet emission. The compact object 
(a neutron star or a black hole) accretes matter from a compan- 
ion star, which can be a massive early type star or a low-mass 
late type star. Since the first discoveries, one decade ago, micro- 
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quasars have raised a strong astrophysical interest because they 
provide small scale counterparts of the highly energetic outflows 
seen in AGNs and quasars. Their short characteristic time-scales 
allow us to study the physics of accretion and outflow in a much 
faster way than in AGNs. Furthermore, because of their binary na- 
ture, one can derive accurate orbital parameters and masses for the 
accreting compact object and its donor star, which may be cor- 
related with jet properties. There are currently around 15 micro- 
quasars known (Ribo 2005) from a populat ion of ~300 X-ray bi- 
naries van Paradiis & van den HeuvellExxj 1200 ll) . but only 
a few with well determined system parameters. Detailed reviews 
on microquas ars can be found in Mirabel & Rodriguez 1 1999|) and 
lFendeJi200l . 

LS 5039 is the optical counterpart of the X-ray source 
RX J 1826. 2— 1450 and it was proposed a s a High Mass X- 
ray Binary (HMXB) by iMotch et all ll997t). Based on archival 
NVSS data and their own VLA observations lMartf. Paredes & Ribol 
( 1998) identified its radio counterpart, which appeared to be a 
persistent non-ther mal radio sourc e. Subsequent VLBA observa- 
tions, performed bv lParedes et alJ feoOOl) . revealed relativistic ra- 
dio jets at milliarcsecond scales, which qualified LS 5039 as a mi- 
croquasar. In the same work, Paredes et al. draw a possible con- 
nection with the unidentified gamma-ray source 3EG J 1824— 1514, 
and suggest, for the first time, that microquasars could be sources 
of h igh-energy 7-r ay emission above 100 MeV (see lParedesll2005l 
and iRibo. Combi & Mirabel 120051 for reviews on other associa- 
tions). Moreover, LS 5039 has been associated very recently with 
the very high energy gamma-ray source HESS J1826— 148, de- 
tected at energies above 25 GeV, rein forcing the association with 
the EGRET source 1 Aharonian et al. 2005). The optical compo- 
nent of the LS 5039 system has been c lassified as a 06.5 V((f)) by 
IClarketalj|200ll) and McSwain et al. <200lLl2004 hereafter M01 
and M04, respectively) derived the first orbital parameters, with a 
proposed orbital period of 4.4267 days and a highly eccentric orbit 
of e = 0.48 ± 0.06. In addition, evide nce for intrinsic pola rization 
at the 3 percent level was presented in Combi et al. (2004) and in- 
terpreted as Thomson scattering in the stellar envelope. The present 
paper reports novel spectroscopic observations of LS 5039 that sup- 
port a revised orbital period of 3.9060 days, a higher mass function 
and significantly lower eccentricity than previous claims. Our re- 
sults have important implications for the nature of the compact ob- 
ject, the presence of CNO products in the optical companion, the 
mass loss of the supernova (S N) explosion, the r ecoil velocity and 
runaway nature of LS 5039 (see lRibo et all2002l and M 04), as well 
as for the i nterpretation of the existing X-ray data (see Reig et al. 
l2003llBosch-Ramon et alll2005l and refere nces the rein) and very 
high energy gamma-ray data lAharonian et all2005|) . 

The observational details are outlined in section|2| Section l3~T1 
presents the analysis of the radial velocities and orbital parame- 
ters while section l3~2l focuses on the rotational velocity calculation 
of the optical star. The determination of the stellar parameters and 
chemical abundances are described in section l3~3l In section|4]we 
present the implications of our new orbital solution for the nature 
of the compact object. Based on all the updated parameters we con- 
duct a general discussion in section|5| 



2 OBSERVATIONS AND DATA REDUCTION 

We observed LS 5039 using the Intermediate Dispersion Spectro- 
graph (IDS) attached to the 2.5-m Isaac Newton Telescope (INT) 
at the Observatorio del Roque de Los Muchachos on the nights 



Table 1. Log of the observations. 



Date 


Object 


No. 


Wav. Range 


Disp. 






sped. 


AA 


(A/pix) 


23-26 Jul 02 


LS 5039 


4 


4045^1910 


0.85 


27-30 Jul 02 




67 


3806-5865 


0.63 


30 Jul 02 


HD 168075 B 


1 






31 Jul 02 


LS 5039 


24 










3 


3920-4666 


0.23 




HD 168075 B 


1 








LS 5039 


2 


4497-5224 


0.22 




HD 168075 B 


1 






1-10 Jul 03 


LS 5039 


101 


3806-5865 


0.63 


4-7 Jul 03 




14 


5350-6860 


0.63 



of 23-31 July 2002 and 1-10 July 2003. A total of 196 spectra 
were obtained with the combination of the 235 mm camera and 
the R900V grating which provided a useful wavelength coverage 
(free from vignetting) of AA3900-5500. We used integration times 
of 300-900 s. The seeing was variable (1-2 arcsec) during our 
runs and we used a 1.2-arcsec slit which resulted in a resolution 
of 83 km s -1 (FWHM). In addition, we also obtained five spectra 
of LS 5039 with the holographic grating H2400B with the aim of 
measuring the rotational broadening (v sin i) of the optical com- 
panion's absorption lines. The spectral resolution of these spectra 
was 30 km s _1 and we varied the central wavelength in order to 
fully cover the range AA3900-5100, where most of the prominent 
He I and He II lines lie. The standard star HD 168075 B, of spectral 
type 07 V((f)), was also observed with the same instrumental con- 
figurations for the purpose of rotational broadening analysis. Fur- 
thermore, we obtained 14 Ha spectra at 58 km s _1 resolution on 
the nights of 4-7 July 2003 to derive the mass-loss rate in the op- 
tical star. These spectra were also taken simultaneously with RXTE 
observations with the aim of studying the correlation between X- 
ray flux and mass-loss rate proposed by Reig et al. 1 2003) and con- 
firmed by M04. These results are discussed in Bosch-Ramo n et all 
(2005). A full observing log is presented in TableQ 

The images were de-biased and flat-field corrected, and the 
spectra subsequently extracted using conventional optimal extrac- 
tion techni ques in orde r to maximize the signal-to-noise ratio of 
the output 1 Home 1986). CuAr and CuNe comparison lamp images 
were obtained every 15-30 minutes, and the A-pixel scale was de- 
rived through 4/6th-order polynomial fits to 53/86 lines (depending 
on the set-up), resulting in an rms scatter < 0.04 A. The calibration 
curves were interpolated linearly in time. 



3 ANALYSIS 

3.1 Radial Velocities and Period Search 

All the spectra were rectified, by fitting a low-order spline to 
the continuum, and re-binned into a uniform velocity scale of 
83 km s . We show in Fig.Qthe averaged spectrum of LS 5039 
along with the template star HD 168075 B of spectral type 
07V((f)). The ratio He I A4471/He II A4541 is somewhat steeper 
for LS 5039, which suggests a slightly earlier spectral type than 
Q 7 V, in good agreem ent with the 06.5 V spectral type proposed 
bv lClarketail fcOOD. Therefore, a synthetic 06.5 V template was 
generated for the cross-correlation an alysis. This was comp uted 
using the NLTE library OSTAR2002 iLanz & Hubenvll2003l) for 
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Figure 1. Normalized average spectrum of LS 5039 (top, and offset 0.4 
for display purposes) and the 07V((f)) star HD 168075 B (bottom). Main 
Balmer and He lines are indicated. Note the steeper He I A4471/He II A4541 
ratio for LS 5039, indicative of a slightly earlier spectral type than 07 V. 
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Figure 2. Radial velocities of the 196 R900V spectra of LS 5039 obtained 
when all the Balmer and He lines are included in the cross-correlation. A 
4-day modulation is clearly depicted. 



T cft = 39 000 K, log g = 3.85, usini = 113 km s" 1 (see 
Sects. l3~2l and 13.31 and degraded to our instrument resolution of 
83 km s _1 . Every R900V spectrum of LS 5039 was then cross- 
correlated with the synthetic spectrum, after subtracting the contin- 
uum. The main IS absorption lines and bands at A3934 (Ca II K), 
A4430, A4501, A4726, A4762, A4885, A5449 and AA5487-5550 
were masked, leaving all Balmer and He lines for the analysis. We 
present the resulting velocities in Fig. [2] Note the night-to-night 
variability which supports a ~4-day periodicity, in line with the re- 
sults of M04. Some nights also show evidence for superimposed 
short time-scale (~hr) variability that is likely produced by a con- 
taminating wind component. 

In order to obtain the orbital period of LS 5039 we have per- 
formed a timing analysis on a database including our radial veloc- 
ities (196 data points) and the ones reported in M01 and M04 (54 
data points), to get advantage of the long baseline and to break the 
1-year alias of our window function. We have noted a small system- 
atic offset of —6 km s _1 between McSwain's velocities and ours 
(obtained by computing the running mean of the two databases) 
and this has been corrected before the analysis. Because the or- 
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Figure 3. Period search results obtained after combining our INT radial ve- 
locities with the ones from M01 and M04. Top: PDM spectrum showing the 
deepest minimum at 0.256 cycle d - 1 . Middle: Scargle power spectrum with 
the strongest peak also at 0.256 cycle d — 1 . Bottom: Periodogram obtained 
with the CLEAN algorithm, showing compatible results. 



bital modulation may not be sinusoidal due to a possible eccen- 
tricity (see M01 and M04), we have employed the phase disper- 
sion minimization (PDM) algorithm I Stellingwernll9 78j). that is 
better suited for non-sinusoidal periodicities (see, e.g.. lOtazu et alJ 
l2002l l 20041) . For compl eteness we hav e also used the standard Scar- 
gle Fourier Transform lScarglel ll982) and the CLEAN algorithm 
(Robe rts. Lehar & Dreher| |l987). All the periodograms were com- 
puted in the frequency range v — 0.05-3 cycle d with reso- 
lution 1 x 10~ 5 cycle d _1 . We show in Fig. [5] the power spec- 
tra obtained by the three methods and they all provide significant 
peaks at a frequency of 0.2560±0.0001 cycle d" 1 , corresponding 
to 3.906±0.001 d. We plot in Fig.|4|all the radial velocities folded 
on the 4.4267 d period of M04 and our favoured 3.906 d period. We 
note that although our radial velocity measurements taken on July 
2002 were compatible with the ephemeris of M04, as it is stated 
in their paper, the new measurements acquired on July 2003 are 
not, therefore ruling out the 4.4267 d period. There are several INT 
data points that deviate from the general trend of the radial velocity 
curve with the 3.906 d period around phase 0.45, which correspond 
to the excursion around HJD 2 452 827.5 visible in Fig.|2] The spu- 
rious data point of McSwain around phase 0.4 corresponds to the 
first value reported in M04, which was obtained from a single spec- 
trum in that night and could be the result of a short-term excursion 
as the one quoted above. Since the radial velocity curve is not si- 
nusoidal, we subsequently decided to fit an eccentric orbit model 
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Table 2. Orbital solutions from an eccentric fit to different data sets. The All solution stands for a fit performed on our INT plus M01 and M04 data for all 
lines. The Balmer, He I, and He II solutions have been obtained with a fixed P or b and To on our INT data alone. 



Parameter All Balmer He I He II (adopted) 

P orb (days) 3.90603 ± 0.00017 3.90603 (fixed) 3.90603 (fixed) 3.90603 (fixed) 

T (HJD-2451000) 943.09 ±0.10 943.09 (fixed) 943.09 (fixed) 943.09 (fixed) 

e 0.31 ±0.04 0.34 ±0.04 0.24 ± 0.06 0.35 ± 0.04 

w(°) 226 ± 8 226.7 ± 3.3 228.0 ± 4.2 225.8 ± 3.3 

7 0ms- 1 ) 8.1 ±0.5 9.2 ±0.8 12.9 ± 0.8 17.2 ± 0.7 

Jfifkmr 1 ) 19.4 ± 0.9 30.4 ±1.6 21.1 ± 1.3 25.2 ±1.4 

aisini(R ) 1.42 ±0.07 2.20 ±0.12 1.58 ±0.10 1.82 ± 0.10 

f(M) (Mq) 0.0025 ± 0.0004 0.0094 ± 0.0015 0.0035 ± 0.0007 0.0053 ± 0.0009 

rmsoffit(kms- 1 ) 6.8 11.2 11.5 9.1 



3.906 d 




0.5 1 1.5 2 

Orbital Phase 

Figure 4. INT (filled circles) and McSwain's (open circles) velocities folded 
on the 3.906 d and 4.4267 d periods using T =HJD 2451 943.09. Two or- 
bital cycles are shown for better display. As can be seen, the new data rule 
out the M04 period. 

(see lWilson & Devinnevlll97ll and Ivan Hamme & Wilsonll2003h 
to our database. As McSwain's velocities have no error bars we as- 
signed equal weight to all the velocities in the fit. Our best solution 
yieldsP orb = 3.90603±0.00017 d and T =HJD 2451 943.09±0.10 
(i.e. time of periastron passage). The rest of the orbital elements are 
listed in the first column of Table|2| 

1 The original Wilson & Devinney code has suffered major up- 
grades since its first release, including significant improvements in 
the underlying physical models. The most recent version of the 
code together with the relevant documentation can be found in 

|f tp ://astro.ufl. edu/pub/wilson /lcdc2003| 
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Orbital Phase (P=3. 90603 d) 



Figure 5. INT velocities folded on the 3.90603 d period with 
To=HJD 2451943.09. The radial velocities have been obtained through 
cross-correlating all Balmer (top), He I (middle) and He II lines (bottom). 
We superimpose the best eccentric orbital fits, the parameters of which are 
quoted in Tablel2l 

Since hydroge n lines in early-t ype stars may be contaminated 
by wind emission JPuls et alJll9961) . we decided to treat different 
line species separately. Therefore, we extracted sets of radial ve- 
locities for the H I, He I and He II lines from our INT spectra. Fig- 
ure[3]presents the radial velocity curves for the three groups of lines 
folded on the 3.90603 d period and T =HJD 2 451 943.09. We have 
also fitted these velocity curves with eccentric orbital models, but 
using the data errors for the weighting scheme and fixing P or b and 
To to the above values. Table|5|lists the best-fitting parameters for 
the separated groups of lines, and the corresponding orbital solu- 
tions can be seen overplotted in Fig. [5] We note a trend in the 7- 
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Figure 6. Relative orbit as seen from above of the compact object around 
the optical component, which lies in the ellipse focus, computed using our 
adopted solution from Table l2l The coordinates are in units of the orbital 
semi-major axis. Relevant phases such as the periastron, apastron, and con- 
junctions are indicated, with the dashed line joining periastron and apas- 
tron. The grey circle represents the optical companion to scale in the case 
of Ro = 9.3 Rq, while the long-dashed ellipse represents the radius of the 
Roche lobe at each orbital phase for the adopted masses of Afrj = 22.9 Mq 
and Mx = 3.7 Mq (see Sect.|4j. The star is at 85 per cent of filling the 
Roche lobe during the periastron passage. 

velocities towards redder values as one moves from Balmer to the 
higher excitation He I and Hell lines. Blue-shifted velocities are 
usually considered a signature of stellar wind, due to P-Cygni con- 
tamination of the photospheric profiles, and can be conspicuous in 
the low members of the Balmer series and some He I lines, such 
as A4471. Furthermore, the fitted solutions to the Balmer and He I 
lines present significantly larger rms scatter than the He II solution. 
This is also consistent with a tentative scenario of variable wind 
contamination in Balmer and He I lines. Therefore, we decided to 
give more credit to the He II solution, which will be adopted as 
the true orbital elements of LS 5039 for the remainder of the pa- 
per. As a test, we have also computed a combined solution using 
the weighted mean of every orbital parameter and we find this to 
be fully consistent with the He II solution within errors, except for 
the systemic velocity which is obviously on the blue side. Our K- 
velocity yields a larger f(M) and compact object's mass range than 
in previous works (M01, M04). To better illustrate the resulting or- 
bit, we show in Fig. [6] the relative motion of the compact object 
around the optical companion as seen from above (i.e., for an ob- 
server with i = 0°). 

3.2 Rotational Broadening 

As a first appro ach, we have fo llowed the technique applied to 
1A 0620-00 bv lMarsh. Robinso n & Wood| < 19941) and described 
in their paper. Essentially, we subtract different broadened versions 
of the 07 V((f)) template HD 168075 B from our Doppler corrected 
average spectrum of LS 5039 and perform a % 2 test on tne resid- 
uals. The 07 V((f)) spectrum was broadened by convolution with 



the rotational profile of iGravl 1 19921) which assumes a linearized 
limb darkening coefficient e. We have taken e = 0.23 which is ap- 
propriate for the stellar parameters of our star (T e « = 39 000 K, 
log g = 3.85, see Sect. I3.3i in the B-band. We performed the 
analysis independently for the two groups of spectra at different 
resolutions, 83 and 30 km s _1 , and find minimum \ 2 f° r broad- 
enings of 105 ± 2 and 126 ± 12 km s _1 , respective ly. Since th e 
template also has a rotational velocity of 79 km s _1 IPennvl 19961) . 
we need to sum this quadratically to get the true rotational velocity 
of LS 5039. This yields v sin i = 131 ± 2 and 149 ± 12 km s _1 for 
the two resolutions. These numbers can be compared with previ- 
ous determinations of 131 ± 6 km s _1 (M01) and 140 ± 8 km s -1 
(M04). 

Alternatively, we attempted a refined deter mination following 
the Fourie r technique describe d bv lGravl ll992l) : for a recent appli- 
cation see iRover et alj 120021) . This technique takes advantage of 
the fact that the Fourier transform of the rotational profile has ze- 
roes at regular positions that depend on the projected rotational ve- 
locity. The whole line profile, being a convolution of different con- 
tributing profiles in the wavelength domain, conserves these zeroes 
in the Fourier domain. While the natural profile does not introduce 
additional zeroes, other broadening mechanisms may do, although 
at the large rotational velocities of the O stars the firsts zeroes of the 
transformed profile are expected to be due to rotation. However, as 
the Stark effect strongly dominates the line profiles of Balmer and 
He II lines, we have only used the He I lines for the v sin i deter- 
mination. This method yields v sin i — 113 ± 8 km s _1 , with the 
error obtained from the dispersion of the individual lines. This is 
a significantly lower projected rotational velocity than any of the 
previously obtained ones. However, this lower value is consistent 
with the finding that a part of the line broadening usually attributed 
to rotation is probably due to some kind of turbulence (see for ex- 
ample Rva ns et all2002t) and, therefore, we give more credit to this 
latter determination. We note that for the fit to the line profiles (next 
section) we have used a projected rotational velocity of 150 km s~ 
in good agreement with M04. 

3.3 Stellar Parameters 

We have determined the stellar parameters using the latest 
version of FASTWIND JSantolava-ReVj Puis & Herrerol 1 19971 
iRepolust, Puis & Herrerdl2004tlPuls et all2005l) . This is a NLTE, 
spherical, mass-losing model atmosphere code particularly opti- 
mized for the analysis of massive OBA stars. We have fitted si- 
multaneously the following lines: Ha, H/3, H7, H<5, He I A4387, 
A4471, A4922, and Hell A4200, A4541, A4686. From the analysis 
we obtain T cff = 39 000 ± 1 000 K and log g = 3.85 ± 0.10. 
While our temperature is consistent with that of M04 our gravity 
is slightly lower than the one adopted by these authors (although 
adopted errors allow for a marginal agreement). Note, however, 
that their gravity, which was determined by model fitting only to 
the wings of the H7 profile, is likely biased upward by the blend 
of N III lines (see details in M04). We also confirm that LS 5039 
is contaminated by CNO products (see M04), with strong N and 
weak C lines, as compared to stars of similar spectral type. It can 
therefore be classified as a member of the ON type group stars, i.e., 
its complete spectral type is ON6.5 V((f)), as given by M04. 

Adopting the best-fitting extinction parameters given by M04 
(i.e. E(B -V) = 1.28 ± 0.02 and R = 3.18 ± 0.07 which leads 
to Av = 4.07± 0.11, different from what they quote), the updated 
My = —4.77 ± 0.15 for an Q6.5V sta r from the new calibra- 
tion by |MartinsISchaerer¥Him^ and V = 11.33 ± 0.02 
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and V = 11.32 ± 0.01 mag from lClark eta?] fcOOll) we derive a 
distance of 2.54 ± 0.04 kpc. We have also considered the slightly 
different values of V reported in the literature, ranging fr om 11.20 
to 11 . 39 mag with typical la uncertainties of 0.03 m ag l Drilling 
ll99lHLahulla & Hiltonll993lMarti et alll998l:lMarti et all2004l) 
and computed the corresponding distances. The weighted mean and 
standard deviation of all these values provides a more realistic dis- 
tance estimate of d = 2.5 ±0.1 kpc, that will be adopted for 
the remainder of the paper. With these numbers we get a ra dius of 
Ro = 9.3±£:6 R © ( see - e.g. lHerrero. Puis & Naiarrcl2 002). From 
this radius we derive a luminosity of log(Lo /L©) = 5.26 ± 0.06 



and a mass of Mq = 22. 9Z 



Mr. 



The mass-loss rate has been derived from the Ha p rofile, 
which appe ars to be variable on a se cular time-scale iReig et all 
120031: M04: lBosch-Ramon et aiT2 005). Our 2003 data show a very 
stable equivalent width with a mean value of EW = 2.8 ±0.1 A 
and no variations above a 5 per cent level over 1 orbital cycle. 
Despite this, we have selected two extreme profiles showing min- 
imum and maximum EW for the calculation, that we have per- 
forme d by using again the late st version of FASTW IND l Pu is et all 
l2005t see examples of use in ReDolust et al. 200^). For the lower 
state (i.e., largest absorption) we obtain a wind mass-loss rate of 
3.7 x 10 -7 Mq yr _1 . While we can put an upper limit for this state 
of 5.0 x 10~ 7 Mq yr _1 , the lower limit is much more uncertain, 
as the sensitivity of Ha to mass-loss rate falls off rapidly at such 
low values. For the high state, we obtain a best value of 7.5 x 10 -7 
Mq yr _1 , i.e., a factor of two if we interpret the small difference 
in the profiles as due to the stellar wind and not to the low signal- 
to-noise ratio. Again, 5.0 x 10~ 7 Mq yr _1 can be considered a 
lower limit for this state, while we can set an absolute upper limit 
to the high state mass-loss rate of 1.0 x 10 -6 Mq yr" 1 . With these 
values, the Modified Wind Momentum of the up per state coincides 
with the values quoted by Repol ust et alj j2004), being a factor of 
two lower for the low state, but still within their uncertainties. We 
also note that our values are a factor ~4-8 higher than those de- 
rived by M04 from data acquired in 2002, but compa re well w ith 
average mass-loss rates in 06.5 V stars (Howarth & Prinia 1989). 



4 MASS OF THE COMPACT OBJECT 

Our revised mass function and stellar parameters provide new con- 
straints on the mass of the compact object in LS 5039. First of all, 
considering the adopted value of /(M) = 0.0053 ± 0.0009 M 
and using an inclination of i — 90° we obtain a strict lower limit 
for the mass of the compact object as a function of Mo- This limit 
is represented by the lowest solid line in Fig. [7] where we also plot 
similar relationships for other inclination angles. For the range of 
possible masses of the optical companion, 22.9^2 9 we obtain 
a lower limit on the mass of the compact object in the range 1.34— 
1.61 Mq, with an additional linear la uncertainty of 0.11 Mq due 
to the mass function error bar. 

Another constraint comes from the fact that no X-ray eclipses 
are seen in LS 5039. With the new e phemeris reported here, 
the BeppoSAX observations performed bv lReig et alj 1200*31) cover 
phases between 0.969 and 0.205, while an eventual X-ray eclipse 
should be centered at phase 0.058 (see Fig. [(J. In fact, periastron 
takes p lace at t = (1.0 ± 0.9) x 10 4 s in figure 1 of Re ig et alj 
(2003), while superior conjunction of the compact object takes 
place at t — (3.0 ± 0.9) x 10 4 s, where the errors come from 
the uncertainty in our To. The X-ray flux is nearly constant dur- 
ing this last interval, implying that X-ray eclipses can be defini- 



7.0 



6.0 



5.0 



4.0 



3.0 



2.0 



1.0 



0.0 



- 






- 


7 


15°/ ! 








/ yA 




-•-"30" 








4 i----" 
















90" \ 











5 10 15 20 25 30 35 40 45 50 

M o ( M sun) 

Figure 7. Mass constraints for the two stars in LS 5039 derived from our or- 
bital solution. The thin solid lines have been calculated by using our adopted 
mass function, f(M) = 0.0053 Mq, and the quoted inclination angles. 
The dashed line indicates the lower limit of Mx from the absence of X-ray 
eclipses. The dotted lines enclose the area of valid solutions for the interval 
of likely values of Mq . The thicker line represents the valid solutions by 
assuming pseudo-synchronization, while the dark-grey region represents all 
their possible Icr errors, and the light-grey region their 3<r errors. The cross 
indicates the case of Mo = 22.9 Mq, which implies ip_s = 24.9° and 
M x = 3.7 M . 



tively ruled out in LS 5039. This condition has been used to com- 
pute the lower limit on AIx as a function of Mo (which in turn 
is also a function of Ro). This is represented as a dashed line in 
Fig.0 which yields i < 64.6° and Mx > 1.49 M for the case 
of Mo = 20.0 M , and i < 63.3° and M x > 1.81 M© for 
Mo = 26.3 Mq. Therefore, this condition provides a compact ob- 
ject mass above 1.49 ± 0.11 Mq (the errorbar due to the /(M) un- 
certainty). This v alue is still below the 1 .75-2.44 Mq neutron star 
mass m VelaX-1 !Ouai ntrellet"aiT2 003). or the (95 per cent confi- 
dence) lower limit of 1.68 M for at least one of the two pulsars in 
binari es with eccent ric orbits within the globular cluster Terzan 5 
iRansom et alj20o"5T) . Therefore, all the previously obtained values 
are compatible with a massive or even a canonical neutron star. 

A further constraint can be obtained by assuming that the op- 
tical companion star is pseudo-synchronized, i.e. its rotational and 
orbital angular velocities are synchro nized at periastron passage 
lKopalll978llClaret & Gimenezll993T) . We have estimated the syn- 
chronization time -scale of the o ptical component of LS 5039 using 
the fo rmalism by IZahnl ll989h (as formulated in IClaret & Cunhal 
1997) and the tidal evolution parameters of the stellar models of 
IClarelH2004l) . We find a synchronization time-scale of about 1 Myr. 
This value is similar to the (rough) upper limit to the age of LS 5039 
(Ribo et al. 2002) and, therefore, the system appears to have had 
time to reach orbital pseudo-synchronism. Note, however, that both 
the synchronization timescale and the age are only accurate as an 
order of magnitude estimate and, consequently, we can only state 
that both quantities do not exclude each other. On the other hand, 
the theoretical orbital circularization time-scale comes out to be a 
much larger value of 10 Myr, compatible with the observation of an 
eccentric orbit. Therefore, assuming pseudo-synchronization (P-S) 
we can combine our determination of v sin i and radius to estimate 
the binary inclination ip-s, i.e. 
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vsini = 2-7T F Ro sinip_s f or bi 

where F = (1 + e) 1/2 /(i - e) 3/2 = 2.22 ± 0.17, and we find 
values of ip_s in the range 26.8-23.2°, depending on Afo. These 
inclinations set strong constraints to the compact object's mass, as 
is displayed by the thicker solid line in Fig. [7] The possible val- 
ues of A/x range from 3.14 to 4.35 M . The la uncertainty re- 
gion (obtained through propagating errors in ip_s and f(M)) is 
marked by the dark-grey area in Fig.0 This leads to compact ob- 
ject masses in the range 2.75-5.00 M©. The 3a uncertainty region 
is also indicated in Fig.0with a light-grey area and yield masses in 
the range 2. 14-6.78 M . The central value of Mo = 22.9 M Q im- 
plies j P _s = 24.9 ±2.8° (la uncertainty) and M x = 3.7±^;g M 
(la uncertainty in all involved parameters), and is indicated by a 
cross in Fig.0 Therefore, if the plausible assumption of pseudo- 
synchronization is correct, the compact object is consistent with 
being a black hole. 

Since there is no significant Ha emission and the binary sys- 
tem does not exhibit stron g X-ray outb ursts during the periastron 
passage (see, e.g., Bosch-Ramon et al. 2005), the donor star is not 
expected to overflow its Roche lobe at any orbital phase. To check 
this, we have generated a grid of solutions in the following way: 
for each possible optical companion mass in Fig.Qwe have iterated 
for all possible inclination angles and obtained the corresponding 
M x (and its error) through the mass function. Then Mx is used 
to compute the Roche lobe at periastron and its uncertainty (using 
the semi-analytical expression in Eggleton 1983). We find that, for 
the possible range of optical star masses Mo = 20.0-26.3 Mq 
and within la, the Roche lobe is not overflown at periastron for 
A/x < 8M Q (or i > 13.4°). 

A final constraint on the inclination angle comes from the ro- 
tational velocity. The break-up speed, or critical rotational veloc- 
ity, can be o btained by means of v CI it = y / (2/3)GM /R (see 
Porter 1996). Considering the limits of our likely values for the 
mass and radius of the optical companion we obtain n cr it = 540- 
580 km s" 1 . Finally, by using our measured rotational velocity of 
v sini = 113 ± 8 km s~ , we can constrain the inclination angle 
to be above 11-12°, which in turn provides upper limits for the 
compact object mass of 8-10 Mq (depending on Mo). 

In summary, the lack of X-ray eclipses, Roche lobe overflow, 
and break-up speed constrain the inclination angle to be in the range 
13-64°, and the compact object mass to 1.5-8 Mq. If the primary 
in LS 5039 is pseudo-synchronized the inclination angle is 24.9 ± 
2.8°, and the compact object mass Mx = 3-7t\[l M Q- 

The inclination angle could be tested through accurate pho- 
tometry near the periastron passage. We have simulated the light 
curve of the tidally distorted optical component using our orbital 
parameters and predict a ~0.02 mag increase in brightness near 
the periastron for i = 45°, which could be revealed by 2-3 mmag 
accuracy photometry. If this variability is not detected or is below 
0.01 mag, the inclination angle must be below 30°, and the mass 
of the compact object above 3.0 M , confirming its black hole 
natur e. We have looked for this effect in our existing photometry 
jMartf et alj|2004l) but the results are not conclusive due to noise 
and scarce phase sampling. More observations are currently under- 
way. 



5 DISCUSSION 

In this section we evaluate the implications of all the updated pa- 
rameters on several properties of LS 5039. 



5.1 Optical companion 

The contamination of the optical companion by CNO products 
may be a consequence of mass transfer of CNO-processed mate- 
rial from the supernova progenitor prior to the explosion (M04). 
Alternatively, it might be caused by mixing processes (induced 
by rotation) which brings CNO processed material from the in- 
ner re gions into the stellar atmosphere jHeger. Langer & Wooslevl 
2000). Indeed, if the system is pseudo-synchronized, then our mea- 
sured v sini = 113 ± 8 km s _1 translates into a large rotational 
velocity of v = 268 ± 34 km s _1 , which is ~0.5fJ C rit • There- 
fore, although we cannot discard the hypothesis of mass transfer 
prior to the SN explosion, we favour the mixing scenario to ex- 
plain the presence of CNO products. We also note that the opti- 
cal companion in LS 5039 is slightly under mass i ve for its spectral 
type, since from the calibration of Martins et al. (2005) an 06.5 V 
star should have A/ spcc = 29.0 Mq, while our measured mass 
of Mo = 22.9^2.9 M © 

covers spectral types between 07 V and 
08.5 V in their calibration. This effect has been observed in other 
HMXB (see, e.g. Jvan den Heuvel et alll983l) . 



5.2 X-ray variability and accretion/ejection energetic balance 

An interesting consequence of our orbital solution is the smaller 
eccentricity, e = 0.35 ± 0.04, compared to McSwain's solution, 
e = 0.48 ± 0.06. A lower eccentricity is easier to reconcile with 
the weak orbital X-ray modulati on, with the X-ray flux varying a 
factor of 2.5, observed by Bosch-Ram on et all 120051) . In a simplis- 
tic Bondi-Hoyle scenario (i.e. spherical accretion through winds) 
our eccentricity, together with the updated masses, provides an X- 
ray variability with a factor of 12 between maximum and minimum 
flux, w hich is still a factor ~5 higher than observed (see lReig et alj 
2003 for details about the method). Even when considering that all 
possible uncertainties conspire in one or the opposite direction, the 
expected variability is in the r ange 7-24, which is a fact or 3-10 
higher than observed. This led Bosch-Ram on et all l2005h to pro- 
pose the existence of an accretion disc that would partly screen the 
compact object from direct wind accretion. 

We can now compare the luminosity released in the vicinity 
of the compact object (i.e., excluding the contribution from the 
companion which dominates in the optical domain) with the avail- 
able accretion luminosity. The first one includes the radio, X-ray, 
high energy, and very high energy (VHE) gamma-ray domains. It 
is maximum in the high energy gamma-ray domain since, using our 
new distance estimate of 2.5 kpc, we have: I/ ra dio, o.i-ioo GHz = 
8 x 10 30 <Martiet alJll998h . L x , 3-30 kev = 1.0-2.5 x 10 34 
iBosch-Ramon^t^lTl200^r after removing t he diffuse backgroun d 
contamination), L 7 , >ioo MoV = 2.7 x 10 35 lHartman et al. 1999), 
and Z/vhe, >250 Gov = 4 x 10 33 erg s _1 lAharonian et all 2005). 
Regarding the accretion luminosity, we have used the Bondi-Hoyle 
scenario to obtain a rough estimate of the accreted matter and the 
energy that can be released from the accretion process. By using 
the updated parameters (Mo = 22.9 M , Ro = 9.3 R , Mx = 
3.7 M , Rx = R Sch = 2GAf x /c 2 = 2.95(M X /M Q ) km, 
P orb = 3.90603 d, e = 0.35, M = 5 X 10~ 7 M Q yr" 1 ), u inf = 
2440 km s~ from M04 and j3 = 0.8 we obtain an accretion lumi- 
nosity averaged through the orbit of < L acc >= 8 x 10 35 erg s _1 . 
Therefore, around 1/3 of the accreted luminosity would be radi- 
ated within the relativistic jet of LS 5039, while the remaining 
available luminosity could be lost by the advection of matter to- 
wards the black hole. Of course, the scenario is not so simple. First 
of all, we should include the kinetic luminosity of the jet, given 
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by Lk = (r — l)Afj ct c 2 (where F is the bulk Lorentz factor of 
the jet, and we have neglected the energy needed to abandon the 
potential well). For a jet with a velo city = v/c = 0.2 (from 

cos 9 = 0.17 ± 0.05 found by Paredes et al. 2002 and assuming 
9 = ip_s ^ 25°), or F = 1.02, we obta in L k = 1.0x l0 36 e rg s _1 
(see details on how to estimate Aij Parede s et alJl2000h . This 
luminosity is slightly higher than the accretion luminosity, although 
this discrepancy could be solved by increasing the mass of the com- 
pact object up to 5 M Q and/or the mass-loss rate of the primary 
up to 10 -6 Mq yr _1 . However, the Bondi-Hoyle accretion sce- 
nario is an oversimplification of the real accretion processes that 
take place in this source, because: 1) the presence of a thick ac- 
cretion disc around the compact object seems to be necess ary to 
launch the relativistic jets in microquasars (see, e.g. Fende3 Eo05l 
and references therein); 2) the presence of a disc also appears to be 
necessary to explain the weak orbital X-ray variability in LS 5039 

1 Bosch-Ra mon et all2 005). Moreover, we note that due to fast ro- 
tation of the optical star and its proximity to fill the Roche lobe 
during periastron passage, we cannot discard additional mass-loss 
in the equatorial plane of the binary system, which could provide 
the needed amount of additional accretion luminosity. We note that 
we have neglected the energy stored in magnetic fields through all 
this discussion. In any case, we can say that the current estimates 
of accretion/ejection luminosities agree quite well by using our new 
set of parameters. We note that if the compact object were a neu- 
tron star with Mx = 1.4 Mq and Rx = 10 km we would have 
< iacc >= 5 x 10 s4 erg s~ 1 . This value is more than one order 
of magnitude lower than for the case of a 3.7 M black hole, and 
~2.5 times smaller than the gamma-ray luminosity of LS 5039. 

5.3 Before and after the SN explosion 

Using our new radial velocity of the binary system 7 = 17.2 ± 

0. 7 km s _1 and an updated proper motions estimate of jj, acoB s = 
4.8 ± 0.8 mas yr~\ fig — —10.9 ±0.9 mas yr _1 (including a new 
radio position obtained with VLA+Pie Town observations; Marti, 
Ribo & Paredes, in preparation), we can recompute the total sys- 
temic velocity of LS 5039. We fin d tWs = 126 ± 9 km s _1 (see 
details on the method in Ribo et al. 2002). On the other hand, our 
improved masses and eccentricity have an impact on the formation 
history of LS 5039. Tidal forces act to circularize the binary orbit 
and hence the current eccentricity e = 0.35 ± 0.04 can be taken as 
a lower limit to e pos t-SN, the post-SN eccentricity. In the context 
of a symmetric SN explosion, e pos t_sN is related to the mass lost 
in the SN event AM through AM = e post _ S N x (Mx + M ) 
which yields AM > 9 ± 2 M and P rc _ c irc = 3.2 ± 0.2 d (us- 
ing our current value s for P or b and e). Wit h these numbers, and 
using the equations in Nele mans et alj ll999h . we obtain a theoret- 
ical recoil velocity of 130 ± 20 km s~ , in good agreement with 
our new space velocity. The previous discrepancy between these 
two values reported by M04 vanishes thanks, mainly, to the lower 
value of e, since the equations are only slightly sensitive to Mx- 
Therefore a high mass loss of ~9 M during the SN explosion 
could provide both the eccentricity and the space velocity that we 
currently observe in LS 5039. We note that our predicted orbital 
period prior to the SN explosion of 1.8 d yields a Roche lobe radius 
of ~ 8.7 Rq, which is compatible within errors to the O star ra- 
dius. Finally, the relatively large mass of the SN progenitor, around 
13 Mp), is c ompatible with the < 10-15 Mq upper limit found by 

iFrver & KalogeralfcOOlli . 

The kinetic energy of the binary system is 4.2± 0.8 x 10 48 erg, 

1. e., merely 4 x 10~ 3 times the energy of a typical SN, although 




Figure 8. VHE gamma-ray fluxes of HESS J1826— 148, the proposed 
TeV counterpart of LS 5039, folded with our orbital ephemeris (P or b = 
3.90603 d with T =HJD 2451943.09). The dashed line represents the 
weighted mean of all data points. Despite the large error bars, there seems to 
be a quasi-sinusoidal variation with a maximum around phase 0.9, a similar 
behaviour to the one seen in X-rays. 

nearly one order of magnitude higher than in GRO J1655— 40 
jMirabel et aHEo 02X On the other hand, from the pulsar birth ve- 
locity value of ~ 400 ± 40 km s" 1 ( Hob bs et alJl2005h . and as- 
suming a pulsar mass of 1.4 Mq, we obtain that their typical linear 
momentum is ~ 560 ± 56 Mq km s . In contrast, the linear mo- 
mentum of the LS 5039 binary system is 3 350 ± 450 Mq km s , 
well above the previous value and above the value of any indi- 
vidual pulsar. To our knowledge, this value is also significantly 
higher than in all other binary systems with measured velocities, 
being a possible exception 4U 1700—37. For this sys tem, moving 
at 70 ± 5 km s" 1 , and usin g different masses from lAnkav et alJ 
1200 ll) and lClarketaill2002h . we obtain a linear momentum in the 
range 2300^-200 Mq km s _1 , encompassing the value found for 
LS 5039. 

We note that once the orbit is re-circularized, with a period 
of 3.2 d, the expected Roche radius will be 15 Rq, well above 
the current radius of the star. After leaving the main sequence, 
the star will fill its Roche lobe, leading to unstable mass trans- 
fer that will prob ably turn-off the X-ray bin ary and microquasar 
phase of LS 5039 iFrank. King & Raind2002h . From the evolution- 
ary tracks bv lMevne^taLfTl994l) we know that the main sequence 
lifetime is between 4.5 and 6.5 Myr (assuming a solar metallicity of 
Z — 0.02 and for stars with initial masses between 40 and 25 Mq). 
On the other hand, the LS 5039 trip from the galactic midp lane to 
its cu rrent position would take ~0.5 Myr (explained in Ri bo et all 
2002 and independent of the distance to the source). Consequently, 
the X-ray binary could survive as a microquasar up to 4—6 Myr 
from now, placing the system at galactic latitudes of —10 to —13°. 
Therefore, relatively nearby systems similar to LS 5039 could be 
the counterparts of unidentified EGRET sources with \ b\ > 5°. 

5.4 Orbital behaviour of the TeV counterpart 

lAharonian et alj l2005h have recently reported the detection of 
a very high energy gamma-ray source at TeV energies, namely 
HESS 11826—148, with a position consistent at the 3a level with 
that of LS 5039. Moreover, the spectral energy distribution at high 
energies makes the association between the HESS and EGRET 
sources virtually certain. These authors state that no periodic varia- 
tions are apparent when folding the data using the orbital ephemeris 
of M04. We show in Fig. [8] the HESS data folded using our new 
orbital period. Despite the large error bars, we notice a possi- 
ble flux variation of a factor ~3 with the orbital period, with 
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a quasi-sinusoidal pattern and maximum around phase 0.9. This 
behaviour is reminiscent to the one recently found in X-r ays in 
RX J1826.2-1450/LS 5039 by Bos ch-Ramon et alJ <2005h . with 
a flux variation of a factor ~2.5 and a maximum around phase 0.8. 
Although further HESS observations would be necessary to con- 
firm this orbital variability, this similarity reinforces the association 
between LS 5039 and HESS J1826-148 and, therefore, practically 
confirms the association between the microquasar and the EGRET 
source. 



6 SUMMARY 

We have reported new optical spectroscopy of the microquasar 
LS 5039 and obtained a new orbital solution. In particular, we find 
Porb = 3.90603 ± 0.00017 d, e = 0.35 ± 0.04, systemic ve- 
locity 7 = 17.2 ± 0.7 km s _1 and a mass function for the com- 
pact object f(M) = 0.0053 ± 0.0009 M Q , significantly different 
from previous results. We have also derived a new distance esti- 
mate of d — 2.5 ±0.1 kpc and a mass of the optical compan- 
ion of Mo = 22.9^2 g Mq. Using this information and assuming 
pseudo-synchronization we obtain an inclination of i = 25 ± 3°, 
which yields to Mx = S.7t\'o Mq. This strongly suggests that 
the compact object in LS 5039 is a black hole. With our new orbital 
parameters there is a good agreement between the accretion and 
ejection luminosities around the compact object. The space veloc- 
ity of the binary system is also in good agreement with the theoret- 
ical recoil velocity in a symmetric SN explosion with a mass loss 
of ~9 M© . Finally, the orbital variability of the TeV counterpart is 
reminiscent to the one seen in X-rays, reinforcing the association 
between LS 5039, HESS J1826— 148 and 3EG J1824-1514. 
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